Calculation of quantum discord for qubit-qudit or N qubits 
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Quantum discord, a kind of quantum correlation, is defined as the difference between quantum 
mutual information and classical correlation in a bipartite system. It has been discussed so far for 
small systems with only a few independent parameters. We extend here to a much broader class 
of states when the second party is of arbitrary dimension d, so long as the first, measured, party 
is a qubit. We present two formulae to calculate quantum discord, the first relating to the original 
cntropic definition and the second to a recently proposed geometric distance measure which leads 
to an analytical formulation. The tracing over the qubit in the entropic calculation is reduced to a 
very simple prescription. And, when the d-dimensional system is a so-called X state, the density 
matrix having non-zero elements only along the diagonal and anti-diagonal so as to appear visually 
like the letter X, the entropic calculation can be carried out analytically. Such states of the full 
bipartite qubit-qudit system may be named "extended X states", whose density matrix is built of 
four block matrices, each visually appearing as an X. The optimization involved in the entropic 
calculation is generally over two parameters, reducing to one for many cases, and avoided altogether 
for an overwhelmingly large set of density matrices as our numerical investigations demonstrate. 
Our results also apply to states of a iV-qubit system, where "extended X states" consist of 



(2 



N+2 



1) states, larger in number than the (2 



JV+l 



I) of X states of N qubits. While these 



are still smaller than the total number (2 — 1) of states of N qubits, the number of parameters 
involved is nevertheless large. In the case of N = 2, they encompass the entire I5-dimensional 
parameter space, that is, the extended X states for N = 2 represent the full qubit-qubit system. 

PACS numbers: 03.65.Ta,03.67.-a 



I. INTRODUCTION 



Quantum correlations in a bipartite system AB have 
enerated interest for various tasks such as computing 
l[, imaging [2[ and metrology Like entanglement, 
which has been proven to be a resource for such tasks 
[l|, quantum discord [|| was proposed to be a resource 
behind the efficiency of the so-called DQC-1 model Q. 
Though this has since been disputed @ , it is of interest to 
understand quantum discord and its potential usefulness 
in quantum technologies. 

Reliable and convenient measures of entanglement, 
such as concurrence Q and positive partial transpose [1] , 
exist only for AB a qubit-qubit or qubit-qutrit system 
0. Quantum discord, on the other hand, is calculated 
very differently, through computation of mutual informa- 
tion as captured by Shannon entropy Q. While the cal- 
culation of eigenvalues and thereby entropy is relatively 
straightforward, a main difficulty in computing quantum 
discord is to calculate the classical correlation, discord 
being defined by default as all that remains of the quan- 
tum mutual information between A and B after subtract- 
ing out the classical correlation between them. This cal- 
culation involves considering all possible measurements 
on one of the subsystems, say A, tracing over A, and 
then computing the entropy which is subtracted from the 
quantum mutual information of AB. This requirement of 
considering all possible measurements so as to account for 
any classical correlation that may exist between A and 
B calls for an extremization procedure that can be te- 
dious, especially with increasing parameters. Only a few 



results are available with, in particular, hardly any an- 
alytical procedures worked out for general systems. For 
qubit-qubit, the entire space being of 15 parameters, re- 
sults [H| are available for a class of 3-parameter states 
that were recently extended by one of us and coworkers 
to 7 [ll|. This class of states have been called X states 
[l2| because they have non-zero entries only along the 
diagonal and anti-diagonal of the 4x4 density matrix. 

X states of two qubits have already found application 
in many studies of entanglement [l2| and discord 
They have the virtue that many calculations can be car- 
ried out analytically which is always an aide to under- 
standing and application. At the same time, they encom- 
pass a variety of states of interest, separable and nonscp- 
arable, classical and non-classical, so that confining our 
studies to them is not unduly restrictive. Recently, the 
algebraic symmetry structure of X states has been re- 
vealed [lH and we have extended the description to X 
states of N qubits [l5[. Again, these X states embrace 
many different states of interest such as W [l6| and GHZ 



17 1 states of three qubits and Dicke states of general N 

In this paper, we show how our previous procedure for 
calculating discord for N = 2 X states, both the tracing 
over the first subsystem A and the extremization, extends 
rather readily to general N; the final step of extremiza- 
tion is at most over two parameters and even that can 
be avoided in many cases by straightforward algebraic 
evaluation at a few values of those parameters. The key 
observation is that even though there is an explosion in 
the number of parameters contained in the density ma- 
trix at larger N, there being (2 N+1 — 1) of them, the X 



2 



structure breaks into calculations of a series of 2 x 2 ma- 
trices whose eigenvalues can be written down in closed 
form easily. Further, so long as A is a qubit, no matter 
what the dimensionality of the second subsystem B, the 
cxtrcmization over the parameters involved in von Neu- 
mann measurements on it, is also limited and standard 
and can be rendered in a "universal" closed form, the 
same that was already encountered in A = 2 ll|. In- 
deed, our results hold even beyond A states, applying 
to any AB of the form qubit-qudit, however large the 
dimension d of B, so long as the density matrix of AB, 
when written as four d x d block matrices, has A char- 
acter for each of these blocks. The density matrix pab 
need not itself have X character. Such states may be 
termed "extended A" , pab having the form 2 (g> dx , the 
rf-dimcnsional part being an X state. The calculation 
of eigenvalues for the reduced density matrix ps after 
tracing over A breaks into that of a series of 2 x 2 matri- 
ces and of the full density matrix pab to at most 4x4 
matrices, all of which can be obtained analytically. 

The outline of the paper is as follows. In Section |TT] 
we present the definition of generalized A-states for N- 
qubits and of extended A-states. In Section HTT1 we will 
outline two formulae, an entropy-based formula and a dis- 
tance formula for calculating quantum discord. We will 
discuss the reduction of the optimization from a two- 
parameter optimization in the general case to a single 
parameter optimization under certain circumstances and 
also point out when an explicit analytical form can be 
written down for the entropic formula for quantum dis- 
cord. These parameters either define the measurement 
basis used to compute the so-called classical correlations 
in the entropic formula or define the classical states C 
of zero quantum discord for the geometric distance mea- 
sure. Finally, we present our conclusions in Section llVl 



II. X STATES OF QUBIT-QUDIT AND 
iV-QUBIT SYSTEMS 

X states were originally described [l2j for qubit-qubit 
systems as those with non-zero density matrix elements 
only along the diagonal and anti-diagonal. The former's 
three real and the latter's two complex parameters form 
a total of seven independent parameters to describe such 
a state. Later, an underlying symmetry structure was 
pointed out, that these states conform to a su(2) x su(2) 
x u(l) sub- algebra of the full su(4) algebra of a qubit- 



qubit system [141 ]. Thus, other 7-paramctcr density ma- 
trices that do not have the visual structure of the letter X 
but share the same algebra belong to the same category. 
This recognition also led to a natural generalization to 
A-qubit systems [l5[. Starting from a single qubit, when 
all states are trivially A and have the su(2) algebra, each 
successive A is generated by doubling and adding a u(l). 
For A = 2, we have su(2) x u(l) x su(2), and upon go- 
ing to A — 3, the symmetry is of su(2) x u(l) x su(2) 
x u(l) x su(2) x u(l) x su(2), a sub- algebra of su(8). 



The sub-algebra deals with a 15-parameter set of states. 
There are (2 N+1 — 1) parameters for general A. 

Similarly, for a qubit-qudit system, there are (4d — 1) 
A states. Given an underlying sub-algebra of states and 
operators, so long as all operators on the system includ- 
ing the Hamiltonian and even possibly those leading to 
dissipation and decoherence lie within that sub-algebra, 
the evolution stays within A character. And, these ac- 
tions can be calculated analytically. For this reason, A 
states have figured in a wide variety of studies in quan- 
tum information. While not the full system of A qubits 
or of qubit-qudit, they are nevertheless a large parameter 
set that is capable of describing a variety of systems and 
phenomena of interest while remaining tractable. 

For purposes of the next section, where we are inter- 
ested in tracing over one qubit, we generalize to an even 
larger category of states. Reduced density matrices ob- 
tained from such bipartite states have the A character, 
while the bipartite states themselves do not. We will 
refer to such states as "extended A states" . In a bipar- 
tite qubit-qudit system, for d odd, these involve (8d — 5) 
and, for d even, (8d— 1) parameters. For A qubits, with 
d = 2 N ~ 1 , there are (2 N+2 — 1) states. These numbers are 
larger than the corresponding ones for X states. While 
still smaller than {Ad 2 — 1) or (2 2N — 1) that represent the 
total systems in general, they are nevertheless substan- 
tial and, for A = 2, actually coincide with the total of 
15 parameters, thus applying to all possible qubit-qubit 
states. For this important system, therefore, our ana- 
lytical result applies to a general density matrix with no 
restriction on its form or its individual elements. 

The non-zero elements of A-qubit A states are pn and 
Pi 2 N +i-i f° r * = A. For extended A states, in 



addition, 



Pi,2 N - l +i 



and 



Pi,2> 



for i = 1. 



,N/2 



and Pi,i-2 N ~ 1 an d Pi i 2 JV +i-(i-2 JV - 1 ) f° r * = A/2+1, . . . , A 
are non-zero. Viewed as four equal blocks, the density 
matrix has each block appearing as an X. An example 
for A = 3 is 
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In the next section, we will use these states to write 
down formulae for quantum discord explicitly. 



III. TWO FORMULAE FOR CALCULATING 
QUANTUM DISCORD 

To establish the first formula, we note that quantum 
discord was originally defined entropically as Q 



Q{p A ») :=I{p A *)-C{p A ») 



(2) 



3 



where the quantum mutual information is defined as 
l{p AB ) := S(p A ) + S{p A ) ~ S(p AB ). Here S(p) is the 
von-Neumann entropy. The second term in the defini- 
tion is the classical correlation defined as 



C(p 



AD 



) = sup {S(p B ) - S(p 

A, 



AB 



(3) 



where = Y,PiS( Pi ). Here 



Pi 



where the measurement operators Ai ~ [/|i)(i|C/^ are 
chosen by transforming the orthogonal projectors H; = 
i = 0, 1 for subsystem A along the computational 
basis kets \i) by a gen eral unitary transformation U , and 
Pi =tr(Ap AB A)Min|. 

For single qubit measurements, the unitary operator 
U can be parametrized in terms of the Pauli matrices 
as U = tl + iy.a, where the four parameters t. y are 



constrained by unitarity: t 



vi + vl 



i 



. The 



three independent parameters can be redefined into the 
unit vector z defined as 

z = {2(-ty 2 +y 1 y 3 ),2(t yi +y 2 y 3 ),t 2 +yl-y 2 1 -yl}. (5) 

This unit vector is relevant for the choice of measure- 
ment directions i = ±z and alternative sets of the param- 
eters t, y may be defined for other orthogonal directions 
[Toj . For any vector V and direction ±z, it is easy to 
establish the vector identity 



Ai(a.V)Ai = ±{z.V)A l . 



(0) 



Furthermore, A; = U\i){i\U^ = (1 ± z.a)/2. Hence, this 
can be rendered in the standard polar decomposition as 



cos 2 (|) 



■ sin ( 



oi<f> 



| sin 6e~ l 
sin 2 (§) 



(7) 



and A_, its parity conjugate, is obtained by the substi- 
tution (0, <p) -> (7T - 0, 7T + (/)). 

All the algebra involved in calculating the post- 
measurement state in Eq. Q can finally be reduced to 
the following simple prescription. Viewing a density ma- 
trix p AB such as in Eq. ([T]) as a block 2x2 matrix of 
(4) four equal blocks, A + p AB A + is calculated by multiplying 
the four blocks, considered in a clockwise manner from 
the top left, by the elements in Eq. ([7]), that is, cos 2 (|), 
isinfle- 1 *, 



sin (If) and isin0e 1 ^, respectively. Tracing 
out subsystem A means adding together the four blocks 
of this post-measurement density matrix. Note that since 
the eigenvectors of this matrix depend only on the unit 
vector z, they are only functions of two parameters, the 
polar angles (0, <fi) of z. In the two-qubit case, these two 
parameters were alternatively understood as arising from 
four parameters k,l,m and n constrained by two equa- 
tions k + I 



11 



1 and to 2 + n A = klm 
Before examining the supremum in Eq. ([3]), we com- 
ment on the choice of extended X-states with respect to 
evaluating quantum discord. Note that there exist two 
separate issues that make the evaluation of quantum dis- 
cord a difficult task. The first of these is the fact that 
since polynomials of degree n > 4 are not generally solv- 
able by radicals [li| , analytical computation of eigenval- 
ues for 2 (g) d systems is not guaranteed for d > 2 (for 
N qubits, d = 2 N ~ 1 ). However, this is possible when 
restricted to extended A-states, the eigenvalue equation 
involving at most a quartic polynomial. 

For N = 3, the post-measurement state corresponding 
to A + up to normalization, namely A + pA + , is given by 



J 



/ p u cos 2 (I) 



P41 COS 2 (f) 

s*p 51 sm(9) 



V s*p 81 sin(0) 





P22COS 2 (|) / 9 2 3COS 2 (|) 

P32COS 2 (I) /933COS 2 (|) 





s> 62 sin(0) s> 63 sin(0) 

s*p 72 sin(0) s*p73 s'm(8) 





pu cos 2 (I) sp 15 sin(0) 





p 44 cos 2 (I) sp 4 5 sin(0) 

s*p 5 4sin(0) P55 sin 2 (§) 





s*p 84 sin(0) p 8 5sin 2 (f) 




sp26 sin(0) sp2i sin(0) 
sp 36 sin(0) sp 37 sin(0) 






P66 sin 2 
P76 sin 2 







P67 sin 2 
P77 sin 2 




sp 18 sin(0) ^ 





sp i8 sin(0) 
2 (I) 



P58 sin 





Psssin 2 (|) / 



(8) 



where s = e _l< ^/2. The post-measurement reduced den- 
sity matrix for B, given by adding the four 4x4 blocks, 
is of course an X state of two qubits whose eigenvalues 
are easily obtained. 

The second issue relates to the optimization involved 
in computing the classical correlations C(p) in Eq. ([3]). In 
evaluating the supremum, we note that the conditional 



entropy is symmetric in i = ±z. It is therefore symmetric 
under parity inversion of z, namely (0, <fi) — > (n — 8, 7r+0). 
Therefore, the search for an optimum can be restricted to 
half the range of one of the angles, either to < < ir/2 
or to < 4> < 7r - The proposal in [ll[ of restricting 
to these extremal values which provides a formula for 
quantum discord for A-qubit A-states that is fully ana- 



4 



lytical and bypasses optimization works when there is 
only one extremum at 9 = n/2. But when there is 
more than one value for the extremum, they occur in 
pairs at possibly other values of the angles and an ac- 
tual optimization is required. Interestingly, as noted in 
[20j| . such an optimization almost always results in the 
solutions 9 = or 6> = n/2 and = or <fi = ir. 
We support this conclusion through our numerical in- 
vestigation of quantum discord for three qubit X-states. 
We sampled 10,000 random density matrices chosen ac- 
cording to the Hilbert-Schmidt measure and constructed 
X-states using p x = J^iEipEj, where Ei are diagonal 
matrices given by E\ = diag(l, 0, 0, 0, 0, 0, 0, 1), E 2 = 
diag(0, 1, 0, 0, 0, 0, 1, 0), E 3 = diag(0, 0, 1, 0, 0, 1, 0, 0) and 
i?4 = diag(0, 0, 0, 1, 1, 0, 0, 0). Discord was then com- 
puted by full numerical optimization and the optimal 
angles corresponding to each randomly generated state 
was plotted on the Bloch sphere. We note that for 99.5% 
of the states chosen this way, the optimal measurement 
corresponds to 9 being at either poles or at the equa- 
tor, as per (Tlj . These optimal measurements arc pre- 
sented in figure (JT|). Furthermore, when either the off- 




FIG. 1. Points correspond to the optimal measurement 
for 10,000 (Hilbert-Schmidt) randomly chosen two-qubit and 
three-qubit X-states obtained by full numerical optimization. 
To accommodate both plots on one Bloch sphere, we exploited 
parity symmetry which requires only half of the sphere, dis- 
playing {8, 0} for two qubit density matrices (thin, black) 
and {9 + tt/2, <f) + tt} for three-qubit density matrices (thick, 
red). For w 99.8% of randomly chosen two-qubit X-states 
and ~ 99.5% of randomly chosen three-qubit X-states, 8 — 
or 9 — n/2 provides the true discord. 

diagonals of the post-measurement density matrix are all 
the same, or when only one of the off-diagonals of the 
post-measurement density matrix is non-zero, it is easy 
to see that (f> does not appear in the eigenvalues of the 
post-measurement density matrix and may be set arbi- 
trarily at (f> — 0. Hence, the optimization reduces to that 



in a single variable 9, which is almost always cither 9 = 
or 9 = n/2. See figure (U). 

Finally, we turn to an alternative formula for quantum 
discord for 2®d systems, a geometric measure defined as 

sua 

^ 2) (p)=min||p- X || 2 , (9) 

where x is a classical state. Such a classical state, in 
general, can be written as 

x = 5>n^®P.P\ (io) 

i=l 

where cIa is the dimensionality of subsystem A and 11^ 
are its projectors. p\ are density matrices describing 
states of subsystem B. For two-qubit systems, classical 
states are defined in terms of one pi, two angles defin- 
ing the projectors and two sets of three parameters 
defining each pi, a total of 9 parameters. 

The 15 parameters that describe any two-qubit state 
on the other hand, can be described in terms of three 
terms, two local Bloch vectors which are three dimen- 
sional and one 3x3 tensor. They are defined via the 
equations Xi = tr(p * cr,- <g> I), yi = tr(p * I ® <Xj) and 
Tij = tr(p * o~i (g> Uj). These three tensors are defined 
for classical states via the equations e = tr((ip\a\4>)) 1 
s + = tr((pipi + P2P2) * <?) and T t] = tx(p * o t (g) Uj). 
Since ||ej| = 1, there is another independent parame- 
ter, namely t = p\ — p 2 (this t is not to be confused 
with the one earlier in the entropic definition). It was 
noted in Q that there are only 9 independent parame- 
ters of x because only three of the 9 elements of T are 
independent, and can be written in terms of the vec- 
tor s_ = tr((pipi — P2P2) * <?), with the relationship 
T = e sT- Hence the three three-dimensional vectors 
e, s± and t constitute the 9 independent parameters that 
define a classical state. This formula for a classical state 
can be substituted into Eq. ^) and optimized with re- 
spect to the three vectors and parameter t, yielding the 
result for a geometric measure of quantum discord for 
two-qubits namely 

D { 2\p) = l(\\x\\ 2 + \\T\\ 2 -k lrL ^). (11) 

Here fc max is the maximum eigenvalue of xx T + TT T . 

Consider classical states of 2 (g) d systems. For 7V-qubit 
states, this corresponds to d = 2 N ~ 1 . These states are 
determined by t = p\ — P2 , two angles defining the pro- 
jectors H; and two sets of (d 2 — 1) parameters that define 
each density matrix. Note again that the local Bloch 
vectors are 3x1-, (d 2 — 1) x 1-dimensional and a correla- 
tion tensor that is 3 x (d 2 — l)-dimensional respectively. 
They are defined via the equations xi = tr(p * G{ <8> /), 
yi = tr(p * I ® Oi) and = tr(p * tr* (8 Oj), there be- 
ing (d 2 — 1) linearly independent operators Oi defining 
an operator basis for the d-dimensional subsystem with 



5 



the additional property tr(OiOj) — Sij. Here 5i t j is the 



Again, it can be seen that the classical states can be 
defined in terms of three vectors and an additional pa- 
rameter. This is done via e = tr ((tp\a\ip}), and s± = 
tr((pipi ±p 2 p 2 ) * 6) with T = e s T . With \\e\\ = 1, 
again another independent parameter is t = p\ — p 2 . The 
three parameters, t and the three-dimensional unit vec- 
tor e, along with the two sets of d 2 — 1 parameters defin- 
ing the two d 2 — 1-dimcnsional vectors s± , define all the 
elements of the three-dimensional vector x, the (d 2 — 1)- 
dimensional vector y and the 3 x (d 2 — l)-dimensional 
tensor T. Thus, having established the Bloch represen- 
tation of the classical state, once again, Eq. ([9]) can be 
optimized with respect to t, eand s±. Writing \\p — x\\ 2 
as 

l|p-xl| 2 = l|p|| 2 + llxl| 2 -2tr(p X ). (13) 

= ^ (1 + I|f||2 + I|y1|2 + I|T||2) + 
+ ^(l + t 2 + ||s + || 2 + ||s-|| 2 ) + 

- -(l + txe + ys + + STs-), (14) 

where ||T|| 2 = tr(T T T). Once again, we have to opti- 
mize Eq. (|13j) over s±, t and e. We first optimize over 
the first two sets of parameters. As in the case of two 
qubits, since the Hessian is positive and non-singular, we 
compute the derivatives and find the minimum by setting 
them to zero. These equations yield 

S^-i(-« +( ,-0. (H, 

— jp = -(-J/ + S+) =0, (16) 

os+ d 

^I|P-X|| 2 1, rpT -* \_ ~* \ f\ 

__ = 5 (-T e + S _) = 0. (17) 

After substituting the solutions to these equations, which 
yield a global minimum, we get \\p — x|| 2 = (H^H 2 + 
j|T|| 2 - e(xx T + TT T )e)/(2d). This is minimized by e 
being the maximum eigenvector of (xx T + TT T ) corre- 
sponding to its maximum eigenvalue fc max . The resulting 
minima again has the same form as Eq. pT|) namely 

D%\p) = ±(\\x\\ 2 + \\T\\ 2 -k max ). (18) 

This completes the proof that Eq. (|lll) represents he cor- 
rect form for calculating geometric measure of quantum 



Kroneckcr symbol and a general density matrix is hence 
written in the Bloch representation as 



(12) 



I 

discord for arbitrary (2 (g> d)-dimensional systems. We 
note that this proof is completely analogous to the proof 
in @. 

For X-states (whose elements have been defined ac- 
cording to [HI]), this formula can be written in a compact 
way. For two-qubit X-states, x = (0,0, d\o) and 

(aoo aoi \ 
<Ho an . (19) 
d n J 

which leads to the eigenvalues of xx T + TT T to be Ao = 
d 2 + d\ x and 

\ ± = \(\a\ 2 ±^W~^) (20) 

where \a\ 2 = ag + agi + a\ + a 2 ^ and a 2 = 2(aooan — 
aoiaio)- So 

W = max(A , A+) (21) 

(2) 

This leads to the formula for D A (p) namely 

D { l\p) = \{\a\ 2 + Ao - max(A , A+)} (22) 
which can be written compactly as 
Df{p) = imin(|a| 2 , \\a\ 2 + A - l -^W^)- (23) 

which generalizes the formula given for states with max- 
imally mixed marginals in jf|. Note that this formula is 
the geometric analogue of the formulas discussed in [l(| , 
@] and [H. Note that Eq. ([23]) involves the full seven pa- 
rameters defining all two-qubit states defined up to local 
unitaries, whereas fl0| and @ involved three parameters 
and [23d involved five parameters. 

IV. CONCLUSIONS 

Constructing analytic formulae for quantum discord is 
hard for two reasons. The first difficulty relates to be- 
ing able to compute eigenvalues of matrices analytically, 
which is related to finding roots of polynomials. The 
second difficulty relates to the optimization over either 
the set of local measurements or over the set of classical 



p = Yd 1 1 + Y Xi(ji ® 1 + Y y^ 1 ® °i + Y Y T v ai °j 



1=1 j=l 
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states, depending on either the entropic or the geomet- 
ric definition of discord. Here we present two formulae. 
The first formula involves a new class of states, namely 
extended Y-states. All steps up to the optimization can 
be carried out analytically, the optimization itself then 
being at most over two parameters (9,4>). For density 
matrices of a certain off-diagonal structure, one param- 
eter <fi drops out, and even the remaining optimization 
in 9 can be avoided for most density matrices to give an 
easily computed, useful formula as shown by numerical 
tests. The second formula established that an existing 
expression for the geometric measure of quantum discord 
for two qubits was valid for (2 ® d)-dimensional systems. 

We note an alternative study of discord for qubit-qudit 
that used local operation and classical communication 
to reduce density matrices to a two-parameter set [24 [. 
Witnesses for discord in such states and a criterion called 
strong positive partial transpose for a non-zero discord 
have also been proposed [25j. Bounds on a quantity 



called measurement-induced nonlocality that is dual to 
quantum discord were given for 2 ® d systems in prjj . 
A related paper [27j relies on [26[ to study arbitrary bi- 
partitions. The entanglement of formation for a certain 
class of 2 <S> d systems were given in [28| by starting with 
tripartite 2 ® 2 <E> d, and tracing over the last to end with 
an X state for the 2 <g) 2 subsystem. Our study differs in 
considering directly a 2 <g> d system, so long as it is of ex- 
tended X type, and evaluating discord (according to two 
different expressions) without any further restriction on 
the parameters involved, numbering 8d — 1 for d even or 
8d — 5 for d odd. Finally, after completion of our paper, 
we became aware of a paper (29j that comes to many sim- 
ilar conclusions to ours for the two definitions of discord 
but, whereas their study is restricted to qubit-qubit, ours 
considers also qubit-qudit and Y-qubits. 
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